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SYSTEM AND METHOD FOR DETERMINING AIRCRAFT TAPELINE 
10 ALTITUDE 

Inventor 
Flash Parlini 

Field of the Invention 

This invention relates generally to aircraft, and, more specifically, to determining 
15 aircraft altitude. 

Background of the Invention 
For aircraft certification and modeling for simulation, the actual height of the aircraft 
during flight must be determined accurately within a threshold value typically of a few feet. 
This is required in order to pass various certifications. One method for calculating aircraft 
20 height uses a barometric loop (baro loop) of inertial reference sensor data. However, the baro 
loop produces some errors that reduce the accuracy when determining actual aircraft altitude. 

Other current methods use height derived from inertial vertical speed with corrections 
for ambient pressure and temperature. However, these methods may not have desired 
accuracy for these applications. This is due to designs of aircraft inertial systems, which bias 
25 their vertical calculations with pressure inputs from the aircraft's sensors. 
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Therefore, there exists an unmet need to more accurately determine aircraft height for 
post-flight testing and aircraft modeling. 

Summary of the Invention 

Embodiments of the present invention provide a method, apparatus, and computer 
5 program product for accurately determining aircraft altitude, impact pressure, and calibrated 
air speed. The determined results may be used for analysis in certification processes, used for 
building flight testing or simulation models that also may used in certification processes, or 
used for other purposes such as data to be used in a flight simulator. 

According to an embodiment of the present invention, altitude information of an 
10 aircraft is determined based on recorded altitude information generated by an inertial 
navigation system (INS) of the aircraft and altitude information generated by a global 
positioning system (GPS) of the aircraft. A static pressure value is generated based on the 
determined altitude information. 

In one aspect of the invention, the altitude information is adjusted based on known 
15 aircraft position defined by a system other than the INS and the GPS. 

In another aspect of the invention, an integration is performed of a temperature 
adjusted vertical velocity value produced by the INS or a double integration is performed of a 
vertical acceleration value produced by the INS, and the result of either integration is 
adjusted according to aircraft pitch, roll, and yaw. A least squares fit is performed between 
20 the INS altitude information and the GPS altitude information. 

In still another aspect of the invention, impact pressure is generated based on the 
generated static pressure and previously recorded pressure information from a pitot static 
system of the aircraft. Calibrated airspeed is generated based on the generated impact 
pressure and aircraft performance data or a simulation model is built based on the calibrated 
25 airspeed 



. 2 - Black Lowe & Graham ^ c 

BOEM-1230AP 701 Fifth Avenue, Suite 4800 

Seattle, Washington 98104 
206.381.3300 • F: 206.381.3301 



Brief Description of the Drawings 

The preferred and alternative embodiments of the present invention are described in 
detail below with reference to the following drawings. 

FIGURE 1 is a block diagram of an exemplary system for performing aircraft height 
determination; 

FIGURES 2 and 3 are flow diagrams of an exemplary process performed by the 
system shown in FIGURE 1; and 

FIGURE 4 is a perspective view of an aircraft performing a test flight. 
Detailed Description of the Invention 

Embodiments of the present invention provide a system and method for accurately 
determining aircraft altitude for use in a simulation model for testing various aspects of an 
aircraft relative to an actual flight test. Referring now to FIGURE 1, a system 20 is illustrated 
for generating aircraft altitude information in accordance with the present invention. The 
system 20 includes a computer 24 that receives data obtained by components of an aircraft 
30. The aircraft 30 includes a data acquisition system 34 that stores data in a memory 36. The 
stored data is received from an Inertial Navigation System (INS) 40, a Global Position 
System (GPS) 42, a Radar Altimeter (RAD Alt) 44, and various aircraft sensors 46, such as 
without limitation a pitot static system, and aircraft pitch, roll, and yaw sensors. After the 
aircraft 30 has completed a designated test, the computer 24 retrieves the data stored in the 
memory 36. The computer 24 retrieves the data from the memory 36 by a direct connection 
or a wireless connection. In another embodiment, the memory 36 includes a removable 
memory device that includes the stored data. The computer 24 receives the memory device in 
a receiving port. 

The computer 24 is a general purpose computer, such as without limitation a personal 
computer, a laptop, a mainframe, or a hand-held computer. The computer 24 includes 
memory, a processor, various user interfaces, such as without limitation a keyboard, a mouse, 
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and a display. The computer 24 determines aircraft altitude or pressure that the aircraft 30 is 
experiencing at various points in time during a test scenario according to an exemplary 
process described in more detail below with respect to FIGURES 2 and 3. 

Referring now to FIGURE 2, an exemplary process 80 is illustrated for generating 
5 highly accurate aircraft altitude information, pressure information, and calibrated airspeed 
(CAS) for use in post-flight test analysis. The process 80 begins at a block 82 where 
barometric pressure is measured at the ground in the vicinity of where the flight test is taking 
place. The barometric pressure measurement is performed shortly before or soon after the 
flight test, or both. This measurement is taken in order to get an accurate measurement of 

10 barometric pressure throughout the period of the flight test. At a block 84, during the test 
flight, the total pressure is measured by the aircraft 30 at a pre-defined sampling rate. Total 
pressure P T is measured by the aircraft's pitot system. At a block 86, a history of static 
pressure values at altitude during the test flight are suitably generated after the test flight has 
occurred. Generation of a history of static pressure values is described in more detail below 

15 with regards to FIGURE 3. At a block 88, for each sample period of time the static pressure 
Ps is subtracted from the total pressure P T in order to get impact pressure. In one 
embodiment, CAS is determined from the impact pressure using Bernoulli's equation. At a 
block 92, a simulation model is built using the CAS and data from other aircraft sensors 46. 

Referring now to FIGURE 3, generation of static pressure values P s is shown as 

20 performed at the block 86 in FIGURE 2. The exemplary process 86 begins at a block 104 
wherein the computer 24 determines change in height values (Ah) using INS information. 
The change in height Ah is determined by double integrating an INS vertical acceleration 
value or taking a single integration of a temperature adjusted INS vertical velocity value. 
Equation 1 below is an example equation for generating the temperature adjusted INS 

25 vertical velocity value. 
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VZIC = VZI * 504.7446 * TAMB (1) 
[145442.2 -HP] 

where: 

VZI is raw inertial vertical speed in ft/sec; 

HP is pressure altitude in feet; 

TAMB is ambient air temperature in deg. K; and 

VZIC is vertical speed from the INS 40 that has been corrected to give tapeline 
vertical velocity. 

At a block 106, GPS altitude information along a test flight is determined. At a 
block 108, the GPS information is compared to the INS information during an appropriate 
test period. At a block 110, improved height values of the aircraft 30 are determined based on 
the comparison. At a block 112, the determined improved height values are adjusted based on 
a known aircraft location point. At a block 114, static pressure P s is generated based on the 
determined height values. 

The following calculations are suitably performed to correct INS height data (i.e., 
generating improved height values of the aircraft 30). A relatively stable period of flight is 
selected for generating the improved height values. Equations 2 and 3 below are two 
different methods for calculating an INS Ah according to aircraft position. 

DZ2 = (trapezoidal) integration of VZIC from [fit time start] to [now] - 
sin(YAW)*cos(PITCH)*DX - (sin(YAW)*sin(ROLL)*sin(PITCH) + 
cos(YAW)*cos(ROLL))*DY + (sin(YAW)*cos(ROLL)*sin(PITCH) - 
cos(YAW)*sin(ROLL))*DZ (2) 

DZ3 = second order (trapezoidal) integration of AZ from [fit time start] to 
[now] - sin(YAW)*cos(PITCH)*DX - (sin(YAW)*sin(ROLL)*sin(PITCH) + 
cos(YAW)*cos(ROLL))*DY + (sin(YAW)*cos(ROLL)*sin(PITCH) - 
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cos(YAW)*sin(ROLL))*DZ 



(3) 



One method for comparing the GPS information to the INS information is shown in 
Equations 4 and 5. 



where: 

DX, DY, and DZ are the distances between the inertial sensor and the vehicle 
reference points in an appropriate body axes system; 

AZ is vertical acceleration (from the INS 40); 

PITCH, ROLL, and YAW are Euler attitude angles of the aircraft 30; and 

ZPDGPS is the height given by the differential GPS (DGPS) system 42, 
corrected for pitch, roll, and yaw to a reference point. 

A second order least squares fit of DZI versus time is calculated: 

DZIFIT = COZ + C1Z * TFIT + C2Z * TFIT * TFIT 
where: 

TFIT is time. 

In one embodiment, a second order least squares fit equation is used if 4 or more 
DGPS points are available, a first order fit equation is used if 3 points are available, and a 
zero order fit equation is used otherwise. If no DGPS data exists, then COZ=C1Z=C2Z=0. 



DZI = ZPDGPS -DZ2 



(4) 



or 



DZI = ZPDGPS - DZ3 



(5) 



T = [time now] - [start time of inertial vertical fit] 



ZPINTU = DZ2 + COZ + C1Z * T + C2Z * T * T 



or 



ZPINTU = DZ3 + COZ + C1Z * T + C2Z * T * T 
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where: 



T is a running time used in the correction; it can extend beyond either or both 



ends of the fit time period; and 



ZPINTU is the unsynchronized height change produced by this method from 



the INS 40 data. 

A point (in time) for synchronizing to a reference height is identified. In this 
embodiment, the reference height is known terrain, where radar altimeter data is used. 
A synchronization constant at time Tl is calculated: 

DZ1 =ZPREF1 - ZPINTU 
Then, for all times of interest: 
ZPBSfTl = ZPINTU + DZ1 
where: 

Tl is the time of synchronization, Tl is used only to identify the point for the 
above equation above; 

ZPREF1 is the reference height at time Tl, obtained by various means, such 
as without limitation from GPS information, corrected radar altimeter data, 
laser altimeter data, stable ground location, position fix with an known object 
(visual or photographic); and 

ZPINT1 is the DGPS-corrected inertial sensor height of the aircraft 30. 

This synchronization can be repeated for different segments of time in a flight test. 
These segments, along with time segments of height data computed by other means, can be 
combined to give a history of aircraft height for the duration of the test period. 

Static pressure Ps is determined according to Equation 6: 



dfPg) = - Ps* (0.010413) 
d(Z) TAMB 



(6) 
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where: 

Z is the height determined above. Z is the tapeline altitude. 
Referring now to FIGURE 4, a perspective view of an aircraft 200 on approach to 
landing on a runway 202 is shown. This is an example test flight that may be analyzed by 
5 embodiments of the present invention in order to produce accurate altitude values in a 
post-test analysis of the aircraft 200 after having flown a flight path 204 to touchdown on the 
runway 202. A calculation of the height data performed by the process of the block 86 
(FIGURE 3) is performed twice. The first time is at the block 108 (FIGURE 3) when the 
GPS information is compared to the INS information. The initial starting location of the INS 

10 information is irrelevant to actual aircraft height, so therefore it is adjusted to the start point 
of the GPS information. The second adjustment occurs at the block 112 (FIGURE 3), where 
a known aircraft location point is used to further adjust the height values that are determined 
in the block 110 (FIGURE 3). FIGURE 4 illustrates one example for determining a known 
aircraft location point. In the test flight shown, the aircraft 200 is at some threshold distance 

15 from the runway 202 (not to be confused with the runway threshold) where the ground under 
the flight test profile (i.e., flight path 204) has been mapped so that its true surface height is 
known. The computer 24 (FIGURE 1) uses stored radar altimeter information at the threshold 
point to determine the height of the aircraft 200 at the threshold point. The determined height 
information at the threshold point is used in the second adjustment. This process may be 

20 likewise performed for a takeoff from the runway 202. 

While the preferred embodiment of the invention has been illustrated and described, 
as noted above, many changes can be made without departing from the spirit and scope of the 
invention. Accordingly, the scope of the invention is not limited by the disclosure of the 
preferred embodiment. Instead, the invention should be determined entirely by reference to 

25 the claims that follow. 
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